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Abstract 

We consider a one-dimensional lattice model with the nearest-neighbor in- 
teraction V\ and the next-nearest neighbor interaction V2 with filling factor 
1/2 at zero temperature. The particles are assumed to be spinless fermions or 
hard-core bosons. Using very simple assumptions we are able to predict the 
basic structure of the insulator-metal phase diagram for this model. Com- 
putations of the flux sensitivity support the main features of the proposed 
diagram and show that the system maintains metallic properties at arbitrar- 
ily large values of V\ and V2 along the line V\ — 2V2 = 7J, where J is the 
hopping amplitude, and 7 ~ 1.2. We think that close to this line the system 
is a "weak" metal in a sense that the flux sensitivity decreases with the size 
of the system not exponentially but as 1/L a with a > 1. 
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The interest to the theory of one- dimensional systems is only partially related to the 
study of organic conductors and other quasi- ID compounds. Another source of interest in 
the ID physics comes from the variety of problems which are either exactly soluble^ or more 
amenable to computational approach. Their solutions give guidance to intuition which can 
be applied to problems in higher dimensions. 

We consider a ID system on a lattice with the following Hamiltonian: 

H = J^2(a]a j+1 + h.c.) + ^ v \i-j\ n i n j (!) 

We study only the filling factor v = 1/2. In the case of the Coulomb potential V\%-j\ = j\ 
one should maintain neutrality and change n« — *> Tij — v. 

We consider the spinless fermion system at T = 0. One can show that for an odd 
number of electrons N the Hamiltonian coincides with that for hard-core bosons. For even 
iV the fermion-boson transformation requires the change of periodic boundary conditions 
into antiperiodic. The particle-hole symmetry can be shown to require that for even iV at 
v = 1/2 the states with total quasimomenta P, ir — P, —P, and P — tt are degenerate. 

The system under study undergoes structural and insulator-metal (IM) phase transitions 
when the hopping amplitude J is varied. The general point of view is that at small J the 
ground state has a crystalline order and is insulating. In the free-fermion limit of large J 
the system does not have long-range order and is metallic. 

In the case of nearest-neighbor interaction and only then the problem is exactly 
solubleM In this case the structural transition occurs simultaneously with the IM 
transition.! In principle, two separate transitions are not forbidden. Nevertheless, in the 
qualitative arguments below we assume that these transitions are connected to each other 
and occur at the same J c . 

We concentrate here on the IM transition in a model with the nearest-neighbor and 
the next-nearest neighbor interactions, the so-called [Vi, V^-model. It has been studied^ in 
connection with the spin version of the Hamiltonian Eq. ([I]). The IM phase diagram for this 
model has been studied recently in Ref. [| 
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We detect the IM transition by analyzing flux sensitivity!'! 5E = \E P — E a \, where E p 
and E a are the ground-state energies for periodic and antiperiodic boundary conditions. For 
simplicity, we take E a to be the lowest-energy state with the same quasimomentum P as 
E p . 

Starting from the ordered phase at J = and using perturbation theory with respect to 
J, one can show that SE ~ J N at small J and hence falls off exponentially with the system 
size L = 2N. For free fermions 5E = ttJ/L. Thus, the dependence of the product L5E on 
L and J is a nice criterion for detection of the IM transition. We obtain this dependence by 
exact diagonalization technique. 

The idea we want to check here is that the IM transition is closely related to the point 
defect with the lowest energy in the crystalline phase. At finite J the point defect forms 
a band. The transition occurs at such J that the lowest edge of the band comes close to 
the energy of the ground state.! At this point the ground state becomes a strong mixture of 
the crystalline and defect states. This mechanism reminds the idea of zero-point defectons 
proposed by Andreev and Lifshitz.0 

Such a simple picture of the transition implies that the critical value of J is determined 
by the energy Ed of the defect at J = 0. The empirical rule we propose is J c = ftEd, where (3 
is some number. For the exactly soluble problem with nearest-neighbor interaction (3 = 0.5. 
For the Coulomb problem E^ = 2 In 2 — 1 = 0.386. Our computations0 show that for the 
Coulomb interaction J c is between 0.17 and 0.3, which gives 0.44 < (3 < 0.77. In the 2D 
case we have found@ that /3 is approximately in the same interval. 

Using the empirical relation 

J c = 0.5E d (2) 

we can construct the IM phase diagram for the [Vi, V^j-model (see Fig. 1). Note that the 
explicit value of (3 is not important for the qualitative results. We choose (3 = 0.5 to get the 
correct value of J c for the case Vi = 0, where it is known exactly. We show below that this 
is a right choice in a wide range of V\ and Vi- 



Two competing crystalline structures exist in the [Vi, V^-model at J = 0. The structure 
1 is •o«o, where • stands for an occupied and o stands for an empty site. The structure 2 

is MOO. 

Dotted lines in Fig. 1 indicate three regions. At J = the structure 1 has the lowest 
energy in the region I, where A = 2V 2 — V% < 0. The lowest-energy defect in this structure 
has energy —A and represents a shift of an electron to the nearest site. The structure 2 is 
stable in the regions II and III, where A > 0. In the region II the lowest defect has energy 
A and is also a shift of one electron. In the region III another defect wins, which has energy 
V 2 . This defect is a "domain boundary", when a portion of a crystal is shifted one site to 
the right or to the left. Such shift, in fact, produces two domain boundaries simultaneously. 

Eq. (H) gives the dependence J c (Vi,V 2 ) that is shown in Fig. 1 with solid lines. These 
lines separate insulating and metallic phases. To obtain J c (Vi, V 2 ) one should substitute into 
Eq. (0) the proper expression for the minimum defect energy Ed(Vi, V 2 ) at J = in each of 
three regions as discussed above. The lower solid line shows the IM transition associated with 
the crystalline structure 1. The upper solid line shows the same transition for the structure 
2. It consists of two straight lines in two different regions, II and III, which correspond to 
the different types of defects. 

Fig. 2 shows the results of numerical computation of L5E/J as a function of J at fixed V\ 
and V 2 for a system of 14 electrons. The data for smaller sizes are not shown. However, they 
have been used to find the critical value J c by extrapolation to 1/L — > 0. At (Vi, V 2 ) equal 
to (1,0), (0,1), and (1,1) our criterion predicts the transition at J c = 0.5; at (4,1) it predicts 
J c = 1. These values are indicated by the points a, b, c, and d in Fig. 1, and by arrows in 
Fig. 2. The value J c = 0.5 is exact for the point (1, 0).i~i The results of extrapolation give 
predicted values for the first three points with a 15% accuracy^ For the point (4,1) we have 
gotten J c = 1.2 ± 0.1. Thus, we may conclude that the Eq. (fj) works very well in a wide 
range of V\ and V 2 - 

The most important prediction of the phase diagram Fig. 1 is existence of a metallic 
region between the solid lines which extends infinitely for arbitrarily large V\ and V 2 close to 
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the line A = 2V2 — V1 = 0. Consider the curves in Fig. 2 corresponding to (Vi, V2) = (1, 0.48) 
and (1,0.52). Now with changing J we are moving almost along the line A = in Fig. 1. 
In the first case we deviate a little towards the Crystal 1, and in the second case — towards 
the Crystal 2. Both lines intersect the IM phase lines at large Vi, V2, predicting J c = 0.02 
in both cases. One can see in Fig. 2 that this prediction is basically fulfilled in the sense 
that the exponential dependence on J disappears near this point. For J > J c the system, 
however, does not look like an ordinary metal, where L5E should be size independent. In 
fact, we have observed a weak dependence of L5E on L in a wide range of J between J = J c 
and J ps 0.4. 

Fig. 2 also shows 5E for (Vi,^) = (1,0.50). Now with decreasing J we are moving 
exactly along the line A = 0. In this case the exponential transition to the dielectric phase 
is absent for arbitrarily small J, in agreement with our phase diagram Fig. 1. However, 
there is some size dependence of L5E along the line A = in the region J <C 1. It can 
be described as SE ~ 1/L a with a > 1. Thus, it is not a regular ID metal where a — 1. 
An alternative interpretation of the same data would be an exponential size dependence 
5E oc exp(— L/£) with anomalously large correlation length £. 

Now we study more carefully the close vicinity of the line A = far from the origin. In 
the region A <C Vi,V2 the spectrum of energies at J = has two scales. The large scale 
is determined by Vi and V2, while the second scale is |A|, which is the energy necessary to 
produce a defect. When A = J = the ground state is macroscopically degenerate. 

To separate these two scales we consider a limit Vi, Vi — > 00, J and A being finite. In 
this limit the size of the Hilbert space can be greatly reduced. Only the states which are 
degenerate at A = J = should be taken into account. These states are such that neither 
three electrons nor three holes occupy adjacent sites. 

The reduction of the Hilbert space size is from C\j 2 to approximately /l-2> where /„ 
denote the Fibonacci numbers, defined by f n = / n _i + f n -2, fo — fi — 1- At large n one 
has@ f n ~ ((1 + V5)/2) n+1 /V5. 

With this reduction we can increase L up to 40 (/3s = .63 x 10 8 ). Fig. 3 shows L5E/J as 
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a function of A/2J obtained for different L. The maximum occurs not at A = 0, as can be 
expected from naive consideration, but at A/2J w —0.6. Accurate size extrapolation shown 
in Fig. 4 demonstrate that at this point 5EL/J stays finite as L goes to infinity. Thus, the 
system at A « 1.2J is a normal metal. The flux sensitivity in the limit L — > oo is less than 
the value 7r for free fermions and is equal L5E/J m 2.5. In the phase diagram Fig. 1 the 
"magic" metallic line A = 1.2J is shown with dashed line. This line appears, obviously, as 
a result of quantum mixture of the two different ordered phases. 

Fig. 3 shows also the energy per particle as a function of A/2 J obtained in the same 
limit. We have not found any singularity in the energy in the region of interest. The gap 
between the ground and the lowest excited states with the same total quasimomentum at 
the magic metallic line scales to zero linearly in 1/L, as shown in the inset to Fig. 3. Note 
that usually a crystalline phase on the lattice has a finite gap. 

The inset in Fig. 4 shows the reciprocal correlation length l/£ = —dln(L5E)/dL as a 
function of A/2 J as obtained from the slopes of the curves in Fig. 4 at largest L. Note that 
the condition £ < L corresponds to l/£ > 0.25. Thus, we have a real exponential behavior 
for —3 < A/2J < 2. At large negative values of A/2J the ground state of the system is the 
crystal with the structure 2 with a small admixture of defects which are fragments of the 
structure 1. At large and positive A/2J one has the opposite picture. In the intermediate 
region the ground state is a mixture of these two structures. If we extrapolate l/£ in each of 
the exponential regions, we find that it turns into zero approximately at the boundaries of 
the metallic strip, shown by two parallel solid lines in Fig. 1. This is natural, since the naive 
picture which leads to Fig. 1 does not take into account mixing of two crystalline structures. 

The small value of £ in the intermediate region suggests that the size dependence of L5E 
is not exponential near the magic line. This would imply the existence of another phase, 
which may be named a "weak metal." If such phase exists, there should be a phase lines 
which separate the weak metal from the normal metal, where L5E is size independent. The 
inset in Fig. 1 shows schematically the region of the normal metallic phase. This diagram is 
similar to the one obtained in Ref. |6], except it predicts an infinite metallic line in the plane 
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(v h v 2 ). 

Finally, we have shown that a simple rule Eq. (0) provides a reasonable description of the 
phase diagram of IM transition in the [Vi, V^j-model. We have found an interesting metallic 
phase which exists at any small values of J. The ground state of this phase is a mixture 
of two crystalline phases with moving boundaries. The nature of a small deviation of the 
metallic phase from the line A = is not clear. 

We are grateful to John Worlock for reading the manuscript. We acknowledge support 
of UCSB, subcontract KK3017 of QUEST, and support of the San Diego Supercomputer 
Center. 
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FIGURES 

FIG. 1. Phase diagram of [Vi, ^J-model. Solid lines show the diagram as obtained from Eq. (2). 
The dotted lines separate regions I, II, and III. The point a is known exactly; the points b, c, and 
d are checked by computations. The long-dashed lines in the main figure and in the inset show 
the "magic" metallic line. The short-dashed lines in the inset show schematically the region of the 
normal metallic phase where L5E independent of L. 

FIG. 2. Dependence of flux sensitivity in units of J/L on J for different (Vi, V2) for the system 
with 14 electrons as obtained by exact diagonalization. The arrows show the transition points 
predicted by the phase diagram. The dashed line shows the free fermion result LSE/J = ir. 

FIG. 3. Flux sensitivity L5E/J for different L and the ground-state energy E per particle for 
L = 40 as functions of A/2 J. The energy E is measured from the classical energy of the crystalline 
structure 2. The inset shows the excitation gap along the magic metallic line vs. 1/L. All results 
are obtained by exact diagonalization in the limit Vi, V2 — ► 00 • 

FIG. 4. Size dependence of flux sensitivity for different values of A/2 J in the limit Vi, V2 — > 00. 
The inset shows the slope l/£ as obtained from this size dependence at large L vs. A/2 J. The 
slope l/£ can be considered as the reciprocal correlation length when £ < L ~ 40. 
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